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Localized plasma structures, such as thin current sheets, generally are associated with localized 
magnetic and electric fields. In space plasmas localized electric fields not only play an important 
role for particle dynamics and acceleration but may also have significant consequences on larger 
scales, e.g., through magnetic reconnection. Also, it has been suggested that localized electric 
fields generated in the magnetosphere are directly connected with quasi-steady auroral arcs. In this 
context, we present a two-dimensional model based on Vlasov theory that provides the electric 
potential for a large class of given magnetic held profiles. The model uses an expansion for small 
deviation from gyrotropy and besides quasineutrality it assumes that electrons and ions have the 
same number of particles with their generalized gyrocenter on any given magnetic held line. 
Specializing to one dimension, a detailed discussion concentrates on the electric potential shapes 
(such as “U” or “S” shapes) associated with magnetic dips, bumps, and steps. Then, it is 
investigated how the model responds to quasi-steady evolution of the plasma. Finally, the model 
proves useful in the interpretation of the electric potentials taken from two existing particle 
simulations. © 2012 American Institute of Physics. [http://dx.doi.org/10.1063/L4747162] 


I. INTRODUCTION 

Thin current sheets with a thickness of the order of typi- 
cal ion scales (gyroradius or ion inertial length) or smaller 
play a key role in the structure and dynamics of space plas- 
mas. Ample evidence in support of this fact is available from 
in- situ magnetospheric observations 1-11 and from numerical 
simulations and theory. 12-27 Thin current sheets are thought to 
be relevant for magnetic reconnection 5,7,9,10,12,14,15,18,20,22,24 
and possibly as the magnetospheric cause of thin quasi-steady 
auroral arcs. 21,25 ' It is characteristic of thin plasma structures 
that they contain significant electric fields. They are based on 
the difference between the properties of the relevant particle 
species, particularly, the small electron/ion mass ratio. They 
occur on scales on which the effects of finite gyroradii or 
inertial lengths cannot be ignored. Their full understanding 
requires kinetic considerations. Magnetohydrodynamics 
(MHD), which frequently provides reasonable results on 
larger scales, is widely unable to resolve thin magnetic struc- 
tures, although often MHD indicates the presence of a thin 
current sheet by a singularity of the electric current density. 22 

On scales exceeding the Debye-length, typical space 
plasmas are quasi-neutral with ion and electron densities 
being kept approximately equal by electric forces. Those 
tend to become larger for smaller length scales. As already 
indicated above, the electric signature associated with a 
localized magnetic structure is of particular interest for 
magnetosphere-aurora connection. Indeed, several regions in 
the magnetotail are possibly associated with quasi-steady (or 
“monoenergetic”) auroral arcs. ’ ’ More diffuse arcs 
presumably are associated with magnetospheric Alfven 
waves. Details are given in a recent review. 29 In the present 
paper, the focus is on the former class. Typically, above the 
acceleration region of a quasi-steady arc the electric held is 


directed perpendicular to the magnetic held. The correspond- 
ing electric potentials are often described to be of “U” shape 
or of “S” shape. 36-41 In this paper, we develop a model that 
associates the electric potential prohle directly with the mag- 
netic held and density prohles of the local structure. 

In the investigation of localized structures in plasmas, 
two main areas can be distinguished. The hrst is concerned 
with the question of how they are generated. In this area, 
most studies apply numerical simulations. For example, par- 
ticle simulations have shown that the slow externally driven 
evolution of a wide plasma structure under a variety of cir- 
cumstances leads to the formation of thin embedded current 
sheets. 24,26 We will return to the electric signature of such 
sheets in Section VI. 

The other area concentrates on the local structure of 
quasi-steady plasma conditions, using analytical methods 
where available. They are widely based on exploiting con- 
stants of the motion, ’ ’ particularly when the depend- 
ence on one or more coordinates can be ignored. Choosing 
this method, we assume a steady-state collisionless plasma 
that is translationally invariant with respect to the Cartesian 
y-coordinate. Then the canonical momentum P y and the 
Hamiltonian H given by 

P y = mv y + qA, H = ™ (v x 2 + v z 2 ) + i // (1) 

with (A ~ < ? A ) 2 + ^ 

are constants of the motion. Here, m and q are particle mass 
and charge, v*, v z are x and z components of particle velocity, 
A and <fi are functions of x, z and denote the magnetic flux 
function and the electric potential, respectively. The way in 
which Eq. (1) is written emphasizes that due to the y 
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invariance the motion can be understood as a planar motion 
in a force field described by the potential \ jj. Throughout this 
paper, we assume that the electric and magnetic fields have 
no y components, such that they are given by 


The variable P y was replaced by P, as there is no danger of 
confusion. 

The first condition that we impose on g s is the quasi- 
neutrality condition 


B = VA(x,z) x e y , E = — V</>(x,z). (2) 

Under these conditions 

fs — F s (P y ,H s ) (3) 

is a Vlasov equilibrium distribution function for arbitrary pos- 
itive F s . The subscript s denotes the particle species. We 
choose singly charged ions (g, = e) and electrons ( q e = —e), 
e being the elementary charge. (P y is used as an independent 
variable and therefore does not carry a species label.) 

Various choices of F s have been used to describe the prop- 
erties of thin current sheets. Schindler and Bim 2 and Bim 
et al 28 used local Maxwellians in the v x , v z plane, such that 

F s = g s (P y )exp(-H s /T s ), (4) 

where T s is the temperature (with Boltzmann’s constant 
being absorbed into T s ). It was found that in the range 
T[ > T e , typical potential differences of the order of Tije 
could be achieved for sufficiently different choices of gi(P y ) 
and g e (—Py). For gi(P y ) = g e (—P y ), the potentials typically 
were an order of magnitude smaller. The sensitivity of the 
potential to the choice of g s suggests that it would be desira- 
ble to put the P y dependence of F s on a more physical basis. 
This paper makes a step in that direction. We assume a 
plasma structure in which the magnetic field strength 
remains large enough for justifying a power expansion in 
terms of the ratio of the ion gyroradius and the length scale 
of the structure, truncated after the second order. Adopting 
Eq. (4), we impose two conditions on the choice of g s . The 
first is the familiar condition of quasi-neutrality. The second 
condition (addressed as “P-condition”) requires that there is 
a frame of reference where for every field line the number of 
particles that have their generalized gyrocenters on that field 
line is the same for electrons and ions. (In zeroth order, the 
generalized gyrocenter reduces to the ordinary gyrocenter, 
details are explained in the following section.) One may 
imagine that the plasma has been generated in that way by 
some kind of source process. As we will see, this model has 
a number of interesting properties that seem worthwhile 
exploring. It is favorably tested with earlier simulation 
results (Section VI). 


II. THE MODEL 

Particle distribution functions of the form (4), written 
explicitly, take the form 


f s (v x ,v 2 ,x,z,P) =&(/») exp - 


m s (v x 2 + v z 2 ) 
27; 


( P-q s A{x,z )) q^z) 


( 5 ) 


n e = Hi , (6) 

where n s = §f s dv x dv z dP is the density of species s. The sec- 
ond condition is the P-condition 

u e (A*) = Ui (A*)y (7) 


where 


u s (A*) = 


S(P — q s A * )f s dv x dv z dPdxdz . 


( 8 ) 


This condition is based on a generalization of the gyrocenter 
for gyrotropic motion, which under the present symmetry 
would be located on the field line where Vy = 0 or, using 
Eq. (1), where the condition is P — q s A = 0 holds. We gener- 
alize this notion to include systems that admit deviations 
from gyrotropy. If P — q s A = 0 holds we say that the particle 
has its generalized gyrocenter on the field line with flux 
value A. Thus, u s (A *) is the number density of particles of 
species s that in the present frame of reference have their 
generalized gyrocenter located on the field line A=A*. 
Then the P-condition (7) postulates that this number is the 
same for electrons and ions. 

The only equation left to assure selfconsistency is 
Ampere’s law 


-AA = /i 0 /(A), (9) 

where /(A) is the electric current density, after eliminating </> 
via Eq. (6). 


III. THE EXPANSION 

Here, we introduce the expansion for small gyroradii. 
The smallness parameter is the ratio of the ion gyroradius 
and the local scale length L, represented by s s = \Jm s T s j 
( eBL ), so that formally one expands in powers of y/ffT s . 
Further below, the electrons will be set to the gyrotropic 
limit m e — ► 0. 

We start with density, which in an explicit form can be 
written as 


n,= 


2n 

VmJ, 


:exp 


<?*</> 

' T s 


Gj(P)exp 


1 

2m x T > 


( P-qsAf 


dP, 

( 10 ) 


where g s was replaced by G s = g s T s 3 / 2 /m s which is treated as 
finite for vanishing m s to ensure finite density. The Gaussian 
kernel of the integral localizes the integrand near P = q s A. 
This suggests substituting P = q s A + SP and expanding G s in 
powers of SP , which gives after SP integration 

n i =(27t) 3/2 exp (G s (q s A)+ r ^-G s (q s A)"^j +0(e s 4 ), 

( 11 ) 
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where the prime symbol denotes differentiation with respect 
to A. Since odd powers of s s do not contribute, truncation af- 
ter the second order means that the error is of order e s 4 . The 
quasi-neutrality condition (6) can now be used to express </> 
in terms of A, 

exp ^e(/>(A) Q- + = U(A), (12) 


where 


U(A) 


G e {-eA) +^G e {-eA)" 


(13) 


and the error term was suppressed. We now make use of the 
smallness of the mass ratio mjmi by taking the limit 
m e lnii — ► 0, obtaining 


rriiTi Gi(eA)" 
2e 2 Gi{eA ) 


mTj 
2 e 2 V(A) 


V(A) 


Gj(eA) 
G e (-eA ) 


-T e /T\ " 


G e (-eA)Y‘ /T 
Gi{eA) 


(18) 


As this equation holds for arbitrary A* and there is no danger 
of confusion, we replaced A* by A. After dividing by 
rriiTi /(le 1 ) and defining 


Y = V(A) 


( G e (— eA) \ r<l/r 

V G,(M) ) ’ 


(19) 


Eq. (18) assumes the form 


Gi(eA)" _ Y" 
G,-(eA) ~Y' 


( 20 ) 


U(A) 


G,-(eA) f m,T,Gi(eA)" \ 
G e (—eA) V + 2e 2 G,(eA) )' 


From Eq. (12), we then find 

e<t> = T -Ai\nU, T = T e +Ti. 


(14) 


(15) 


Let us now turn to the P-condition (7). Explicitly, the func- 
tion u s reads 


u s (A*) = 


In 

VrnJ? 


G s (q s A*) 


exp - 


, . .^2 q s <t>{A) 


2 nuT, 


(A-A*) - ■ 


V(A)dA, 

(16) 


where V(A) = J 4 ds/B is the differential flux tube volume, 
ds being the arc length differential on field line A and B the 
magnetic field magnitude. An expansion analogous to that 
applied to density gives 


(2 ) 3 / 2 

u s (A*) = — G i (^ J A*)exp 


qs4>{A*) 


(2nf l2 m s T s ( 

+ e 2 e 2 Gs{ChA 4 eXp 


V(A*) 

q s <P{A*j 


v(A*)y 

(17) 


where now the prime symbol denotes differentiation with 
respect to A*. 

After inserting </> from Eq. (15) into Eq. (17) and 
applying the limit of vanishing electron mass, we obtain 
an explicit expression for the P-condition (7). Analyzing 
that condition, we first find that in the zeroth order it is 
satisfied identically. This is intuitively clear, because to 
lowest order the gyrocenter and particle locations coincide 
so that the P-condition is an automatic consequence of 
quasi-neutrality. The next higher nonvanishing order terms 
are linear in rrii (order g/ 2 ). In that order, the P-condition 
(7) gives 


For any 7, considered an arbitrary function of A, 
Eq. (20) is read as a differential equation for G/(A), which 
has the general solution 


G,(A) = F(A) 


c + Co 
\ 4 0 



( 21 ) 


where c,co are arbitrary constants for any fixed Ao. From 
Eqs. (11), (12), and (19), we find 


G,-(A) = 


n{A)V{A) 
(2? t) 3/2 F(A) ' 


( 22 ) 


Equating (21) and (22) gives 


n 0V0 
Fo 2 


A 


+ Co 


dA 

yW 


n{A)V(A) 

m 2 


(23) 


where the factor ( 2n ) 3 ^ 2 was absorbed in the arbitrary con- 
stants, and then the constant c was expressed by /2o,Eo,7o, 
the values of n,V,Y at Ao- Further, introducing non- 
dimensional quantities h = n/no ) V = V /Vo, and Y 
= Y /To, Eq. (23) assumes the form 


A 

1 + do 

A 0 


dA' 

F(A') 2 



(24) 


where do is an arbitrary constant replacing c 0 . 

Equation (24), read as an integral equation for T, has a 
unique positive solution (satisfying the boundary condition 
at A = Aq) determined by 


(. 


Y = nVex p \ do 


A0 


dA' 

h(A')V(A!) 


(25) 


Using the definitions of U (lowest order) and T given by Eqs. 
(14) and (19), we obtain the electric potential </> from Eq. (15) 
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e(<p - (p 0 ) 


A 

Ti h dA' 

~2 y + 0 J h(A')V(A')’ 

A 0 


(26) 


where oc is an arbitrary constant replacing do and 
0o = 0( A o)- 

This equation holds in the frame of reference where the 
P-condition is satisfied. A transformation to a reference 
frame that moves with velocity rj in the y direction gives in 
that frame 


/ / , x Ti h dA' 

«(*-*■) = -2 + , (A , WA ,) 


-ri(A-A 0 ). (27) 


A 0 


Equation (27) is the desired result for two-dimensional con- 
figurations. It is comforting that the knowledge of the ion 
temperature, density, and flux tube volume suffices to find 
the electric potential as a function of A. 

In the special case of one-dimensional fields with trans- 
lational invariance with respect to x and y, the magnetic flux 
function A depends on z alone, such that the magnetic field 
has an x component only. In the derivation of the expression 
for </>, the spatial integrations in the x, z plane is replaced by 
z integrations with the flux tube volume V being replaced by 
1/B. Accordingly, one finds 


instead of Eq. (28). On the right hand side of Eq. (30), den- 
sity n is replaced by pressure p associated with the x, z plane, 
such that h = p. In the ID cases, one can then take into 
account Eq. (9) (in integrated form) by introducing the pres- 
sure balance 


B 2 

P + — =p c , (31) 

2po 

where p c is a constant. Thus, for one-dimensional fields, <fi 
can be expressed in terms of the magnetic field alone. 


IV. EXAMPLES 


Here, we apply the findings of the previous section to 
simple one-dimensional structures. Where appropriate, we 
use z as the independent variable instead of A; note that A(z) 
is monotonic as B does not vanish (the expansion breaks 
down if B becomes too small). The plasma is assumed to 
reach asymptotically homogeneous states for z = ±oo. We 
use quantities B = B/B 2j p = p/pi, etc., where the quanti- 
ties labeled by the subscript “2” are taken at z — > oc (or at 
the right boundary for finite systems), so that B(oo) = 
1, p( oc) = 1. The pressure balance (31) then gives 


p = i+T( 1 ~b 2 ), p 2 


fi 2 2 

2hqP2 ' 


(32) 


e{(p - (p 0 ) 


A 

T 

yin (nB) + a 
A 


B(A') 

h(A') 


dA' — rj{A — Ao). 


(28) 


Another difference in comparison with the 2D case is 
the additional constant of the motion P x , which in the 
absence of the x component of the vector potential can be 
replaced by v x . Thus, 


F„ = gs(Py)ex P 


nhvJ__hA 

2 T u T 2 J ’ 


Ills, rs 

h s = yv z 2 + i> (29) 


is a Vlasov equilibrium distribution function. It takes into 
account anisotropy by attributing different temperatures to 
the motion in x and in y,z. Carrying out the same expansion 
as above for that case, one recovers the expression (28) with 
the only modification that 7} is to be replaced by Tu. The 
temperature associated with the x motion does not enter ex- 
plicitly. Although T\ s appears in n s and u s , it disappears after 
defining g s = y/fn s / (VT\~ s T 2s )G s , where G s plays the same 
role as in the isotropic case. 

The expressions (27) and (28) are the results of power 
expansions and break down when the formal values of </> 
become too large. This can pose a limitation to the choice of 
the free constants a and rj. For instance, if one considers a 
one-dimensional plasma structure where n and B assume 
nonvanishing constant values for z = ±oo, the constants a 
and f] have to vanish so that we find 

e(0~0o ) = -y ln(«B) = -yln(pB) (30) 


Using Eq. (32) in Eq. (30) and setting </>(oo) = 0, we find for 
(p = ecj)/Ti the expression 

cp(z) = -lln^(z) +T(l -B(z) 2 )£(z)). (33) 

As the first group of examples, we consider magnetic dips 
and bumps, assuming that B(z) is symmetric with respect to 
z = 0 and has a single extremum (at z = 0). In that case, 
c p(z) has either a single extremum or three extrema. The 
number of extrema, as well as the extreme values, are deter- 
mined by the parameters B(0) and jS 2 alone. Thus, for B( 0) 
kept fixed, two otherwise different choices of B(z) provide 
different potentials cp(z), but if /i 2 is fixed also they have the 
same qualitative shape, as defined by the number of extrema 
and by the extreme values. This is illustrated by Figure 1. 

Figure 2 gives the qualitatively different regimes of the 
cp(z) profiles in B(0), space. As the figure indicates, for 
suitable parameters a magnetic dip {B{ 0) < 1) can be associ- 
ated with either an electric dip (or “U-shape”) or an electric 
bump. More complex c p profiles with additional side minima 
also occur. The same applies to magnetic bumps (5(0) > 1). 

Next, we consider magnetic steps with monotonically 
increasing B(z) and again 5(oo) = 1. Here, the parameters 
that determine the qualitative shape of the electric potential 
are B\ =B{— oo) (or for finite systems the value at the left 
boundary) and /? 2 - Correspondingly, the two different mag- 
netic profiles of Figure 3 (with same parameters) have the 
same qualitative potential shape. 

Figure 4 gives the different regions where qualitatively dif- 
ferent c p shapes occur. The figure shows the presence of S 
shaped potentials of different orientation with and without 
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FIG. 1. Two examples for magnetic dips. The solid and dashed curves 
shown in the upper graph correspond to two different choices for B(z), 
picked arbitrarily except that they have the same parameter B(0) (in addition 
to symmetry, single extremum and B( oo) = 1). The lower graph shows the 
corresponding electric potentials cp given by Eq. (33) with f$ 2 fixed also. 



FIG. 2. Parameter space B( 0) vs. fi 2 for magnetic dips and bumps as defined 
in the text. Shown are the regions that correspond to different qualitative 
shapes of the electric potential cp{z) 1 sketched inside the boxes. Same quali- 
tative shapes are indicated by the same color. The region above the black 
curve is unphysical (negative plasma pressure), the green line separates 
magnetic dips (B( 0) < 1) from magnetic bumps (B(0) >1), on the red 
curve <p( 0) changes sign, on the blue curve the two side minima on one side 
merge with the central extremum, and on the magenta line the side minima 
on one side disappear when their cp values reach 0. 




FIG. 3. Two examples for magnetic steps. The solid and dashed curves 
shown in the upper graph correspond to two different choices for B(z), 
picked arbitrarily except that they have the same parameter B\ (in addition to 
monotonic increase with z and B (oo) = 1). The lower graph shows the corre- 
sponding electric potentials cp given by Eq. (33) with a fixed value of p 2 . 


modification involving a minimum. Note that on the red line, (p 
becomes U-shaped although the magnetic profile is S-shaped. 

V. SLOW EVOLUTION 

Steady state models are largely interpreted as snapshots of 
slowly evolving systems, driven by a suitable boundary condi- 
tion, such as addition of magnetic flux. This raises the follow- 
ing question. Suppose our model applies to the electric 
potential of a given initial snapshot of such a slowly evolving 
system. Does it then remain valid during the subsequent 
evolution? We will answer this question for the cases of mag- 
netic dips/bumps and magnetic steps as considered in the previ- 
ous section. Conveniently, we refer to earlier analytical studies 
of quasistatic evolution of one- and two-dimensional systems in 
the gyrotropic limit. 24,26 There, it was demonstrated that in that 
limit, the one-dimensional evolution is governed by simple sim- 
ilarity transformations and can be expressed through a single 
parameter, say b(t), representing the external driving. We 
found that for the quantities that are relevant in the present 
context, the transformations are 

B(z,t) = b(t)B 0 (b(t)z) 
n(z,t) = b(t)no{b{t)z) 
p(z,t) = b{tfpo(b{t)z) 

<KV) = 4> 0 (b(t)z), 

where p is the zz component of the pressure tensor, summed 
over species. The subscript 0 refers to the initial state. 
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FIG. 4. Parameter space B \ vs. /i 2 for 
magnetic steps. Shown are the regions 
that correspond to different qualitative 
shapes of the electric potential <p(z), indi- 
cated in the boxes. Steps where B 
decreases with z (i.e., B\ > 1) are left out, 
because these shapes can be reduced to 
corresponding shapes with B \ < 1 by re- 
versal of the z axis and a renormalization. 


Assuming constant asymptotic states (at z = ±oo), we can 
use Eq. (30) for the model potential, which we write as 

^>m( z ) = -£ ln («( z )5( z ))> (35) 

where, as before, the bar symbol indicates normalization 
with respect to the value at +oo, so that </> m ( oo) = 0. The 
time variable is suppressed. 

Within our approximations and by our assumption that 
the model holds initially, the initial potential is given by 

M z ) = <Pmo( z ) = - T f e \n{n 0 {z)B 0 (z)). (36) 

According to Eq. (34) at some later time t, the potential is 

<t>( z ) = <Po( bz ) = -^ln[n 0 (bz)B 0 (bz)). (37) 

Again using Eq. (34) we find h{z) = h$(bz ), B(z) = Bo 
(bz) such that the model at time t given by Eq. (35) becomes 

</> m ( z ) = -^ln(»( z )£(z)) = b<p(z), (38) 

where we used that T scales as b. 

Thus, 

<t>( z ) =\<t>m( z )- (39) 

So, strictly speaking, the model is not conserved during the 
evolution. However, the violation consists only in the pres- 
ence of the constant factor l/b. The adjusted model potential 

(fra = ( 40 ) 

would provide correct modeling (within the present 

simplifications). 


Thus, our model </> m , if adjusted with a suitable multi- 
plier 1 /b, is relevant for an entire slow evolution if that evo- 
lution contains a single state to which (j) m applies exactly. In 
particular, the qualitative potential shape, within present 
approximations, is correctly represented by the model. This 
also follows from the fact that the parameters B( 0), /? 2 of 
Figure 2 and B i , /? 2 of Figure 4 are conserved under Eq. (34). 
In the following, we will refer to our model as the “relaxed 
model” if adjustment of the electric potential by a constant 
multiplier is admitted. 

VI. APPLICATION TO TWO EARLIER SIMULATIONS 

To gain further insight into the usefulness of our model, 
we applied it to two simulations carried out earlier. 26 Sur- 
prisingly, we found that in both cases Eq. (39) holds approxi- 
mately for suitably adjusted values of the parameter b so that 
the relaxed model proves applicable. 

The results are shown in Figures 5 and 6 for the two 
simulations. Plotted are the simulation magnetic field, the 
simulation potential </>, the model potential </> m , and the 
adjusted model potential (j) a . 

Surprisingly, in both cases the relaxed model curves 
(small-scale-broken lines) show fair agreement with the sim- 
ulation (solid line), while the model (j) m (dashed curves) 
deviates grossly from the simulation potential. The parame- 
ter b was found by a best fit procedure, which gave b = 3.6 
for the first simulation (Fig. 5) and b = 6.0 for the second 
(Fig. 6). 

It is instructive to confirm that the potential shapes of the 
simulations can be identified with the corresponding model 
shapes in Figures 2 and 4, respectively. The magnetic 
step of Figure 5 has the (simulation) parameters B\ =0.06, 
p 2 = 0.10, which is close to the red curve in Fig. 4, on which 
(j) m is exactly U shaped, consistent with the approximate U 
shape of the simulation potential. The approximate U shape 
of the simulation potential in Figure 6 has the parameters 
B( 0) = 0.67, P 2 = 0.05, which in Figure 2 lie in the U shape 
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FIG. 5. The figure provides an application of the present model to an earlier 
particle simulation (Section IV of Ref. 26, interval 10 < z < 40). The upper 
graph shows the magnetic field profile. The lower graph gives the electric 
potential of the simulation (solid curve), the model curve (j) m based on simula- 
tion data (dashed line) and the adjusted model potential (f) a (fine-scale broken 
line). The parameter b = 3.6 was determined by a minimum variance fit. 



FIG. 6. This figure corresponds to Figure 5 for the second simulation 
(Section V of Ref. 26). The parameter b — 6.0 was determined by a mini- 
mum variance fit. 


regime between the blue and green curves. Alternatively, we 
may consider the right flank of the U shaped </> profile in Fig- 
ure 6 as a step with parameters B\ = 0.67, = 0.05, 

which, as expected, yields a step in Figure 4. 

VII. SUMMARY AND DISCUSSION 

This paper addresses properties of electric potentials 
that exist in localized steady plasma structures, specifically 
thin current sheets embedded in uniform or weakly varying 
magnetic fields. A two-dimensional model is presented that 
is based on finite ion gyro-radius effects and expresses the 
potential in terms of ion temperature, density, and magnetic 
flux tube volume (Eq. (27)). The model is valid for a particu- 
lar class of distribution functions given by Eq. (4). Further, it 
uses an expansion for small deviation from gyrotropy and 
besides quasineutrality it assumes that electrons and ions 
have the same number of particles with their generalized 
gyrocenter on any given magnetic field line, the P-condition. 
There is no magnetic field component in the invariant direc- 
tion. The one-dimensional specialization gives the electric 
potential in terms of density and magnetic field magnitude 
(Eq. (28)) and if (via pressure balance) the density is elimi- 
nated, the potential is expressed in terms of the magnetic 
field alone. For homogeneous asymptotic states at large |z|, 
the model assumes a particularly simple form (Eq. (30) or 
Eq. (33)). 

One might wonder why the expressions for the electric 
potential, which is the result of an expansion into ion non- 
gyrotropy, do not involve the ion mass. The reason is that 
the zeroth order of the P-condition is satisfied identically, so 
that the next higher nonvanishing order determines the zer- 
oth order of G t . This property is carried over to the potential 
</>, so that it is determined by zero-order quantities alone. 
The electron temperature drops out, because in the present 
limit of m e /mi — > 0, the electrons do not contribute to finite- 
gyroradius effects. 

The ID model is employed to study the electric potential 
forms of magnetic dips, bumps, and steps (Figures 1-4). 
Each of these elementary magnetic structures can be associ- 
ated with a variety of different potential profiles, their quali- 
tative shape depending only on two parameters. A magnetic 
dip can be associated either with a potential dip (“U” shape) 
or a potential bump, transitional forms can also occur 
(Figure 2). A magnetic step can have an electric potential 
that is either “S” shaped or “U” shaped or of a transitional 
shape (Fig. 4). 

Electric potential shapes play a central role in observed 
quasi-steady auroral arcs. U-shaped potentials above the 
acceleration region corresponding to converging perpendicu- 
lar electric fields are frequently observed to be associated 
with downward acceleration of electrons and upward electri- 
cal currents. 34,39 ^ 1 S-shapes and transitional shapes have 
been observed also. For instance, Marklund et al. 34 presented 
spacecraft data indicating a transition from a U-shaped to an 
S-shaped potential. 

If our model applies to an initial state of a slow evolu- 
tion driven by external action, it remains to be applicable 
during the evolution, however, only in its relaxed form 
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allowing the potential to be adjusted by a constant multiplier. 
This does not affect the qualitative potential shapes, which 
are conserved in the entire evolution sequence. 

Density and magnetic held strength have to be bounded 
away from zero to avoid breakdown of the expansion. The 
expansion can readily be extended to higher orders, should 
this be required for particular applications. Our model does 
not cover plasmas for which the mass ratio is not 

small (e.g., electron-positron plasma), but the method can 
easily be generalized to include such cases. 

When, tentatively, we applied the model to two earlier 
simulations (Figs. 5 and 6), it came as a pleasant surprise 
that both cases were found to be covered by the relaxed 
model. How is this possible, given the fact that our model is 
subject to stringent or arbitrary assumptions such as the 
weak gyrotropy condition and the P-condition? Weak gyro- 
tropy means that £/, the ratio of the ion gyroradius over the 
length scale of the magnetic field, must be small compared 
to 1, but this condition is alleviated by the fact that the error 
is of order g/ l 2 3 4 5 6 7 . In fact, the weak gyrotropy condition is rea- 
sonably satisfied in the cases in question; the maximum error 
amounts to 6% for the first simulation (Fig. 5) and is near 
10 -6 for the second (Fig. 6). Regarding the other assump- 
tions, the situation is less clear. It might be significant that 
the agreement between model and simulations applies only 
to the relaxed version of the model which involves a fitting 
parameter that widens the range of applicability of the 
relaxed model compared with the non-relaxed model sub- 
stantially. For example, the agreement could be explained if 
the simulation state can be understood as a member of a set 
of slowly evolving equilibria with one member being cov- 
ered by the non-relaxed model (see Section V). It seems also 
possible that the actual range of applicability is wider than it 
might seem from our derivation. Investigations beyond the 
present scope seem required to provide a firm answer. 
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